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We present a systematic approach based on Bloch vector’s treatment and the Magnus quantum 
electrodynamical formalism to study qubit manipulation by a train of pulses. These investigations 
include one of the basic processes involved in quantum computation. The concrete calculations 
are performed for tunneling quantum dynamics, multiple resonance and off-resonance excitations of 
qubit driven by Gaussian pulses. In this way, the populations of qubit states due to multiple reso¬ 
nant interactions are investigated for various operational regimes including: single-pulse excitation, 
two-pulse excitation with phase shift between pulse envelopes being controlling parameter and for 
excitation with sequential pulses. In the last case, we demonstrate the formation of quasienergetic 
states and quasienergies of qubit driven by train of identical pulses. In this case the transition 
probability of qubit exhibits aperiodic oscillations, but also becomes periodically regular for definite 
values of the quasienergy. 

PACS numbers: 85.25.-j, 03.65.-w, 03.67.-a 


I. INTRODUCTION 

Interactions of quantum systems with pulse train are of 
great interest for many scientific and technological appli¬ 
cations, such as qubit manipulation or quantum comput¬ 
ing, secure communications and ultra-precise measure¬ 
ments. This direction is similar to the method of nuclear 
magnetic resonance for the radio-frequency control and 
manipulation of spin systems and atoms [ll 1^ . In this 
way, the train of pulses with the name of composite pulses 
has been originally developed in a pioneering work 

The method of control of quantum systems by pulse 
train has also been applied in atomic physics (see, for ex¬ 
ample: 0), in magnetometry with solid state quan¬ 

tum sensors Q and molecular spectroscopy Q- Some 
other applications have been noted in Ref. (8|. 

This direction also involves the tunneling dynamics 
of time-dependently driven nonlinear quantum systems. 
The physics of driven quantum tunneling are associated 
with a wide variety of interesting phenomena and effects 

i- 

On the other hand, decoupling methods based on the 
multiple pulse techniques are very perspective for open 
quantum systems, for decoherence control in quantum 
engineering |10l-ll2l|. In this way, experimental applica¬ 
tions have been obtained for trapped io ns Il3l4l5l| , atomic 
ensembles and spin-based devices |l7l| . 

In this spirit, we emphasize the idea of improving the 
degree of quantum effects in open systems, as well as 
obtaining qualitatively new quantum effects through ap¬ 
plication of the sequence of tailored pulses. This ap¬ 
proach was recently exploited for formation of high de- 
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gree continuous-variable entanglement in the nondegen- 
erate optical parametric oscillator [Tsl - f^ . for investiga¬ 
tion of quantum interference in mesoscopic domain |22j |. 
for production of Fock states and for realizing photon 
blockade in Kerr nonlinear resonator driven by pulse 
train [23 - 1^ . and for demonstration of quantum chaos 
at level of few quanta [13, . 

The control of qubits by pulses is a basic process in¬ 
volved in quantum computation. In this direction su¬ 
perconducting qubits are among the most promising sys¬ 
tems for the realization of quantum computation. Re¬ 
cently, Rabi oscillations, manipulation and control of 
states of single superconducting qubit by resonant mi¬ 
crowave pulses have been extensively studied [29l - l33j| . 

Recently, the dynamics of superconducting qubit 
driven by external field with time-modulated amplitude 
and the phenomenon of Rabi oscillations have been con¬ 
sidered |34l | close to the experimental scheme on the 
frequency-modulated transmon qubit [s^. Supercon¬ 
ducting qubits usually have short coherence time, there¬ 
fore to decrease the time for performing gate operations 
a large-amplitude external fields should be applied. The 
dynamics of a qubit driven by large-amplitude external 
fields in the case of driving around the region of avoided 
level crossing has been also studied (see, for re¬ 

views) . 

In the present paper we present a systematic approach 
for investigation of tunneling quantum dynamics of two- 
level systems, particularly, superconducting qubits inter¬ 
acting with a single pulse as well as with a train of pulses. 
In this way, we analyse the Rabi model on base of both 
the Bloch vector’s treatment and the Magnus QED for¬ 
malism for calculation of time-evolution operator. The 
Bloch vector, three-dimensional unit vector, provides a 
convenient pictorial representation of qubit state [s^ 
while the Magnus formalism [s^, Ii3| leads to expansion 
of the time-evolution operator without well known Dyson 
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time-ordering in QED. We formulate the Magnus formal¬ 
ism in so-called Furry representation that allows to ap¬ 
ply the rotating-wave approximation (RWA) as well as 
to consider naturally the system beyond RWA. In addi¬ 
tion, we demonstrate that this approach is also valid for 
investigation of pulse-driven multi-qubit systems. 

This approach allows to obtain general explicit solu¬ 
tions for single-qubit populations in dependence from 
time-dependent Rabi frequency that can be calculated in 
details on base of the Magnus expansion. Note that this 
result cannot been obtained in the framework of Dyson 
time-dependent perturbation theory and is valid for var¬ 
ious interactions including also the regime beyond rotat¬ 
ing wave approximation. On the other side, using the 
Magnus formalism in the Furry representation allows us 
to formulate an effective truncation procedure for Mag¬ 
nus series. So it will be strongy shown that in RWA the 
operator of time-evolution is only determined by the first 
term of Magnus series and effects beyond RWA can be 
calculated as corrections to the resonant part through the 
high-order Magnus terms. 

It is demonstrated in this paper that the results can 
be also obtained for a wide range of frequencies of driv¬ 
ing field corresponding to multiple resonant excitation 
regimes of qubit on one side and involve interactions with 
both single pulse or the train of pulses with arbitrary en¬ 
velopes on the other side. The concrete calculations will 
be done for the train of pulses with Gaussian envelopes 
and with different phases. The duration of pulses, time 
intervals between them and the relative phases of pulses 
are the free control parameters. Thus, this configura¬ 
tion provides more flexibility for manipulation and con¬ 
trol of qubit states, particularly, it allows us to engineer 
the phase-controllable populations of qubit states. Such 
approach is useful for performing a large number of quan¬ 
tum gates that usually are designed with atom like qubits 
interacting with applied driving pulses with finite dura¬ 
tions. 

Another interesting regime occurs when the pulse train 
consists of identical pulses and hence displays the peri¬ 
odicity in time. We will demonstrate that quasienergies 
and quasienergetic states (QFS) (or so-called the Floquet 
states) for the qubit driven by train of identical pulses are 
realized in multiple resonant regimes. Thus, the regular 
oscillations of the state population is realized if the factor 
of the quasienergy and the period becomes equal to tott, 
where to = 1,2,.... 

Note, that at first, the QFS of the composite sys¬ 
tem consisting of an atom and time-periodic e.m. field 
have been considered first time in the papers (dll - l^ . 
These states provide a classical counterpart to well known 
atomic-dressed states in which the coupling to the 
laser is described by a classical field, whereas the coupling 
to the vacuum must be described in second quantization. 

There have been several experiments on nonlinear and 
quantum optics that have been interpreted in terms of 
quasienergy levels including basic experiments on the 
resonance fluorescence and the probe absorption spec¬ 


troscopy for a two-level atom in a strong laser field. QFSs 
for a two-level atom in the dichromatic field have also 
been studied in a series of papers (see, for example, [dsl - 
Applications of QFS and quasienergies to Josephson 
qubits in driving field have been done in several papers 
[ 33 , Idol - lsijl . including revie w p aper on Landau-Zener- 
Stiickelberg interferometry [^, probe spectroscopy of 
QFS application of the Floquet theory to Cooper 
pair pumping [i^, and observation of Stark effect and 
generalized Bloch-Siegert shift in the experiment with su¬ 
perconducting qubit probed by resonant absorption via 
a cavity The experiments on the Rabi oscillations 
in monochromatically driven Josephson qubits have been 
performed and interpreted on the base of dressed states 


The paper is organized as follows. The single 
driven qubit with time-dependent interaction Hamilto¬ 
nians along the z axis is studied in Sec. |n]in the Bloch 
vector’s treatment and by using the Magnus expansion. 
In this section, the general results for the populations of 
qubit states are reported, particularly, for the regimes of 
multiple resonant excitations including: single-pulse exci¬ 
tation and two-pulse excitation with phase shift between 
pulse envelopes being controlling parameter. Moreover, 
the modifications of this approach for multi-qubit sys¬ 
tems, including details for two-qubit system and off- 
resonance interaction beyond RWA are presented in Sec. 
cn The multiple resonant excitations with sequential 
pulses are considered in Sec. Ell This analysis involve 
also manipulation of qubit by train of identical pulses 
as well as formation of QFS and quasienergies. Finally, 
some conclusions are given in Sec. HV] 


II. PULSED DYNAMICS OF A SINGLE QUBIT 
IN THE FURRY PICTURE AND IN THE BLOCH 
VECTOR’S PATH 


In this section, we derive general expressions for popu¬ 
lations of a single qubit as well as we consider the appli¬ 
cation of the approach for multi-qubits systems, partic¬ 
ularly, for coupled two qubits. However, detailed calcu¬ 
lations will be done for a single driven qubit for various 
multiple resonant excitation regimes including: single¬ 
pulse excitation and two-pulse excitation with phase shift 
between pulse envelopes being controlling parameter. 


A. Time-evolution of qubits on base of the Magnus 
expansion. The probability of excitation and the 
Rabi frequency 


At hrst, we consider the qubit driven by a field of pulses 
with arbitrary envelopes on base of the Hamiltonian 

H{t) = -g{t)dz + AA 


( 1 ) 
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where time-dependent interaction takes place along the 
z axis with the coupling constant 

g{t) = ^eo + A{t)cos{ujt + 0). (2) 

Here eg is the electronic energy difference between the 
ground and excited states of the qubit in the absence 
of the driving field, A is the tunneling amplitude be¬ 
tween the basis states, g{t) includes the interaction be¬ 
tween the external time-dependent field and the two-level 
system, and A(t) is the amplitude of external field, 
ay, Gz are three usual Pauli matrices. The Hamiltonian 
(ED describes various physical systems including quan¬ 
tum two-level systems with nonlinear tunneling mecha¬ 
nism driven by external field @ and also superconducting 
qubit or JJ artificial atoms with A being the Josephson 
coupling energy. Up to now, numerous proposals and 
demonstrations of quantum phenomena in superconduct¬ 
ing qubit have been demonstrated in the fields of atomic 
physics and quantum optics (see, for example review pa¬ 
pers [HIH). 

We start with time-evolution equation of the system 

*^|vk(t)) = i7(t)|vk(t)) (3) 

in the Parry picture |'k(t)) = (7|$(t)). The transforma¬ 
tion leads to the equations 

i^=-g{i)^zU, (4) 

= A(7-id,U|$). (5) 

The operator U is calculated as 

U = exp J g{t')dt' ■ Gz^ (6) 


and time-evolution equation of the system can be rewrit¬ 
ten in the following form 




where 


h = A 




(7) 

( 8 ) 


where <7+ = {gx + iGy)j2, d- = (gx — 'idy)l2 and 


m = f g{t')dt’. (9) 

Jo 

Following a standard perturbation approach time- 
evolution of quantum states is described in the Dyson 
expansion form as 


|$(t)) = T exp 



l<i>( 0 )). 


( 10 ) 


where T is the time ordering operator and |$(0)) is an 
initial state of the system. Considering below high-order 
terms in the perturbation expansion we note that the 
Hamiltonian does not commute at different times due to 
Dyson-type time dependent integrations. This circum¬ 
stance usually makes the analysis of high-order processes 
complicated. On the other hand the unitary evaluation 
operator S{t, 0) can be also presented in the other form: 
so-called Magnus expansion (see, for example [1^), which 
seems to be more preferable for some specific systems 
than the standard Dyson expansion. In this represen¬ 
tation the time-evolution operator is a true exponential 
operator 

S{t,0) = T exp ^—i J h{t')dt''^ = exp(D(t)), (11) 

where 11(1 = 0) =0 and a series Magnus expansion for 
the operator in the exponent takes place 

OO 

= ( 12 ) 

k^l 

For instance, the three terms of the series are given as 

Jo 

Mt) = dt"[h{t'),h{t'% (13) 

^3(0 = ^ J dti J dt2 J dt3(^[h{ti),[h{t2),h{t3)]] 

+ [Hh), [h{t2),h{ti)]]^. 

Thus, the Magnus approach provides the approximate 
exponential representations of the time-evolution opera¬ 
tor of the system. It is not difficult to illustrate a connec¬ 
tion between Magnus series and Dyson time-dependent 
perturbative series. For instance, it is easy to check that 
the second-order Dyson term of perturbation theory cor¬ 
responds to the second-order Magnus expansion. It fol¬ 
lows from the following equality 

dti dt2T (ji{ti)h{t2)^ = dtih(ti) 

+ [ dti [ dt2[h{ti),h{t2)]{lA) 

Jto Jto 

and the structure of Di(l) and D 2 ( 0 - 

A simple description of a two-state system presents 
the state vector as a point on a two-dimensional sphere, 
the Bloch sphere. Therefore, we represent the evolution 
operator as an infinitesimal rotation of the Bloch vector 
in abstract space 

\^{t + dt)) = S{t + dt,t)\^{t)), (15) 
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S{t + dt, t) 


= 1 - 



For the system that initially was in the ground state |1), 
( 71 ( 0 ) = 1, ( 72 ( 0 ) = 0, the population distribution P 2 {t) 
(16) of the state |2) is 


through an angle dO = 2/S.dt about the axis directed along 
the unit vector with the components 

= cos(2A(t)), 

riy = sin(2A(t)), (17) 

Uz = 0 . 

Note, that interpretation of the Bloch vector’s path be¬ 
yond rotating wave approximation for a two-level system 
interacting with a monochromatic field has been done 
recently 

It is easy to realize that in general the time-evolution 
operator in the path of Bloch ball can be represented 
in the simple form for all terms of the series expansion 
(US- Indeed, the all terms of the Magnus expansion are 
only expressed through the Pauli matrices dx, dx, dx 
and hence the following representation of time-evolution 
operator can be derived 


S{t,0) = exp(f7(t)) = exp (^—iG{t)d'^ 

= I cos(|G(<)|) - i{p{t)d) sin(|G(t)|), (18) 

where the vector p{t) = G{t)/\G{t)\. In this way, the 
problem is reduced to calculation of the function |G(t)| 
and the vector p as the series of Magnus expansion. In 
the second-order Magnus expansion we obtain from Eqs. 
0 


Gxit) = A f dt'cos{2A{t')), 

Jo 

Gy{t) = A ( dt's\n(2A(t')), (19) 

Jo 

Gzit) = f dt' f dt"sin[2(A(t") - A(t'))], 

Jo Jo 


and G{t) = |G(t)| = ^ Gl{t) + Gl{t) + Gl{t) . Thus, the 
vector state reads as 


|$(f)) = exp (^-iG(t)d^ |$(0)) (20) 

and hence 

|vl/(t)) = U{t) exp [-iG{t)d) l'h(O)). (21) 

On the whole for the probability amplitudes of qubit 
states 


|vI/(t)) = Gi(<)|l) + G2(t)|2) (22) 


P 2 {t) = {l-pl)sm\G{t)). (24) 

We notice that this result is obtained in the most gen¬ 
eral form for arbitrary time-dependent field amplitude 
A{t) and for a wide range of external field frequencies. 
In the general case, however, there is not compact ex¬ 
pressions for the components of the function G{t) that 
are given as the Magnus series expansions. The first 
term of this expansion m coincides exactly with the 
simple exponent solution of time-evolution operator with 
the Hamiltonian h{t). The every other n-th term r2„(t) 
contains multiple integral of combinations of n — 1 com¬ 
mutators containing n Hamiltonian operators h{t). Nev¬ 
ertheless, the result (HI has some advantages from a 
computational point of view. Particularly, it is easy to 
get the simplest expressions for the population of qubit 
in RWA describing interaction of two-level system and a 
single electromagnetic mode for the case when the fre¬ 
quency of mode is near resonance with the qubit transi¬ 
tion frequency and the coupling constant of interaction 
is weak (for atomic systems, see [sgI - Is^ I. 


B. Multiple resonant excitations 

In this subsection the case of multiple resonant interac¬ 
tion with pulses is considered. The resonance condition 
is formulated using the requirement that the oscillating 
terms in time are vanished in the integrals of Eqs. m- 
Thus, this condition is formulated for the frequency uj 
and the electronic energy difference as Eq — Noj = Ajv <C 
Eq. At first, we consider Eqs. (1^ . dH) in the following 
form 

I /■* 

A{t) — -Egt + / A{t') cos{ujt' + 9)dt' 

= h,, + + _ f‘ sinM' + D)3y)^^, 

2 Lj Jq oj at' 

(25) 

for the case of adiabatic pulses for which duration of 
pulses T ^ In this approximation the expression 

(EH) is simplified as 

./.N I / N sin(wt-I-0) 

A(t) =—Eot + A{t) - . (26) 

Z UJ 


we obtain 

Gi(t) = [cos(G(t)) — ipz sin(G(t))]Gi(0) 

- 1 - _ pj^)sin(G(t))G2(0), (23a) 

G2(t) = e~'^^'\-ipx + Py) sm{G{t))Gi{0) 

+ [cos(G(t)) +ipz sin(G(t))]G2(0)(23b) 


Indeed, the last integral in Eq. (1251) can be omitted in 
comparison with the other two terms in the adiabatic 
regime, that it is easy to check for pulses with Gaussian 
envelope 

jt-tkd 


Ak{t) = Age 


(27) 
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Calculating Eqs. (I19|) with the approximated formula 
(051) we also use the well known formulae with the Bessel 
functions. In the result we can obtain 


n 



^ ■) 


(28) 

where Jn{x) is n-th order Bessel function of the first kind. 
For the case of exact 7V-th order resonance, eo = Noj, we 
get 


« J-N 

= (-1)^^7V (29) 

In this approximation the Hamiltonian ([8]) transforms 
to 

hit) = (-l)^A. 

(30) 

Thus, the components of the vector G{t) up to second- 
order of Magnus expansion are calculated as 


Gxit) s 

S (-I)^Acos(A6») 

f Jn\ 

/2A(t')' 

)dt'. 

JQ 

V w , 

f 

Gy it) S 

s (-I)^Asin(iV6») 

f Jn( 

^2Ait')\ 

)dt', (31) 

lo ' 

V / 

Gzit) - 

i 0. 





As we see, these Eqs. dSU don’t involve the terms of 
second-order as has been noted in Introduction. The 
component Gzit) going from the second-order Mag¬ 
nus term and hence equal to zero in RWA approxima¬ 
tion. It is easy to explain this truncation of Magnus 
series in RWA from the general point of view. Re¬ 
ally, the high-order terms of Magnus expansion Ct 2 it), 
^sit), ... involve commutators of Hamiltonian operators 
[h{ti),h{t 2 )], [hiti),[hit 2 ), hits)]], ... that are equal to 
zero for the case of Hamiltonian in RWA (see, formula 
(30)). Thus, in RWA the operator of time-evolution is de¬ 
termined only by lowest exponential term involved hit). 

On the whole for the population of excited state we 
obtain the following result 



FIG. 1. The population of excited state R 2 (t) for the case of 
first-order resonance in dependence of time-intervals in units 
of l/eo- The parameters are: A/eo = 0.3, Teo = 10, tu/eo = 
1. Solid line (a), Ao/eo = 0.19; dashed line (b), Ao/eo = 0.12; 
dotdashed line (c), Ao/eo = 0.25. 



FIG. 2. The population of excited state P 2 (t) for the case 
of second-order resonance in dependence of time-intervals in 
units of l/eo- The parameters are: A/eo = 0.3, Teo = 10, 
(jj/eo = 1. Solid line (a), Ao/eo = 0.35; dashed line (b), 
Ao/eo = 0.25; dotdashed line (c), Ao/eo = 0.5. 


is small in opposite to the first-order resonance config¬ 
uration. As will be demonstrated below the complete 
population inversion takes place for both resonant con¬ 
figurations in cases of multi-pulse excitations. 


C. Qubit excitation beyond RWA: calculations up 
to the second-order of Magnus expansion 



Below we present application of this result for excita¬ 
tion of qubit by a single pulse at the frequency w and 
with Gaussian envelope (see, Eq. (I?71) '). if Ci(0) = 1, 
C2(0) = 0. The results corresponding to the case of first- 
order resonance are depicted in Fig. [T] As we see, the 
probability allows for a region of high transition proba¬ 
bility with increasing the amplitude of single pulse. The 
case of second-order resonance is demonstrated in Fig. 
[2 As we see, in this case the transition probability T 2 (t) 


In this subsection we shortly discuss above intro¬ 
duced approach in the regime beyond RWA, i.e. for the 
Rabi pulsed model in the regimes where the RWA can 
not be applied. In this way, we consider single-qubit 
off-resonance interaction with strong external Gaussian 
pulses. We calculate numerically qubit populations due 
to non-resonant excitation on the general formulas (1231) . 
(p4)) with high-order Magnus terms. As it can be seen 
from Eqs. m, the first-order Magnus terms (propor¬ 
tional to A) are Gxit) and Gy it) components of the vec¬ 
tor Git), while the second-order term (proportional to 
A^) involves only z component, Gzit). On the whole, the 
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frequencies but with different amplitudes, A{t) = Ai{t)- 
A 2 {t). The exponential function in this case reads as 


,,, 1 , ,,sin(a;t + 0i) , , , sinlwt + 6*2) 

Ht) = T;£ot + Mt)— - -^+A2it)^ -^ 

Z (jJ LU 


and the following expression can be obtained 


(33) 






nAijt) 

V ^ 


2^2 (t) 


gi[£o + (iil+’i2)‘^]tg™l®lgin2®2 


■ (34) 


FIG. 3. The population of excited state 1^2 (t) beyond RWA 
in dependence of time-intervals in units of 1/eo for (71(0) = 1, 
^ 2 ( 0 ) = 0 and the parameters: A/eo = 0.45, Aq/eq = 0.4, 
Teo = 3.5. Solid line (a), lj/eq = 0.5; dashed line (b), u/eo = 
1.5; dotdashed line (c), oj/eo = 1.1. 


Magnus expansion up to the second-order involve all com¬ 
ponent of the Bloch vector in Ea. dMl) and third-order and 
high-order terms only are modified these components. As 
shows calculations the small parameter of this expansion 
is the ratio A/eq < 1. Thus, below we study numerically 
the temporal evolution of the qubit state vector (I22D, (ESD 
up to the second-order Magnus terms. Note, that in cal¬ 
culations based on the formulas (ED the terms depending 
on w -h Eo oscillate very rapidly and have been neglected 
in the rotating wave approximation (see, subsection A). 
In this subsection, we will explore the effects of these 
counter-rotating terms beyond RWA on calculation on 
the formulas (IT^l) including these oscillation terms. The 
typical results are depicted in Fig 13] for three cases of 
non-resonance interaction with the frequencies of the ex¬ 
ternal field: w/eq = 0.5, w/eq = 1.1 and uj/sq = 1.5. 
As we see, in this approach beyond RWA the popula¬ 
tion displays time-dependent oscillations due to involving 
counter-rotating terms in qubit dynamics in contrast to 
analogous results for the resonant case depicted in Fig IT] 
The other peculiarity with respect to the resonance case 
is that population of qubit decreases in time during the 
pulse. This effect is stipulated by the negative term pi in 
Fq. (l24l) appearing in the second-order Magnus expan¬ 
sion. The Figl3]also shows an asymmetry in non-resonant 
excitation of qubit state connected with values of the fre¬ 
quency relatively to the energy difference between ground 
and excited states of the qubit. Indeed, the population 
for the case of w/eq = 1.5 exceeds analogous result for 
w/eq = 0.5. It should be mentioned that calculation of 
the third-order Magnus term (see, Fqs. ED) for used 
parameters is only slightly change the obtained results. 


D. Phase-dependent two-pulse resonant excitation 

In this subsection the dynamic of qubit in the result 
of two-pulse excitation is considered. We assume that 
the external field consists of two pulses with the same 


We assume multiple resonant interaction if Eq — Noj = 
0, where ni -I- n 2 = —N. In RWA we reach 

JV 

g2.A(t) ^ jA,(u;(t)) (^^5?) ' 

where w{t) = (Z^ + z'^ — 2zZcos((/>))^/^, z = 

Z = (f> = 9i — 62 + TT. This formula is valid when 

|ze±*‘^f< jZj. 

This formula can be applied for time evolution that 
involves independent interactions with pulses, i.e. there 
is time interval between two interactions subjected by 
pulses, separated by interaction-free interval. Fach of 
the independent interactions are described by a slowly 
varying RWA Hamiltonian (l30t . Nevertheless, we assume 
that pulses depend on the phases of the field envelopes. 
Thus, in this regime the complete expression for expo¬ 
nential function for an arbitrary moment in time reads 
as 


„2iA(t) _ 


= (-l) 


N 


Jn 


2Ai(t) 


UJ 

- Jn 


2 A 2 {t) 

LO 


^-iN 02 


(36) 


On the whole the population of excited state is calculated 
from (IT9t . (I24l) and (l36l) as 

P2(t) = sin^ , (37) 


where 

= Ji(i) + J2 (^) + 2cos(JV( 02 - 9i))ji(t)j2(t) (38) 

and 


jk{t) = Jn dt'. (39) 

This expression takes place for the general case of 
multi-order resonant interactions in term of the factor 
N. As we see, in this configuration the controlling phase 
factor involved in the probability is the phase difference 
between consecutive pulses. 

Now we focus on investigation of phase-dependent ef¬ 
fects for one-photon resonant configuration considering 
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FIG. 4. Transition probability population P 2 {t) for the case of 
two-pulse excitation in dependence of time-intervals in units 
of l/cQ. The parameters are: A/eo = 0.3, Tea = 10, r/T = 6, 
uj/eo = 1. Solid line (a), Aq/eq = 0.19; dashed line (b), 
Aq/eo = 0.12; dotdashed line (c), Aq/eq = 0.25. The phase 
difference of two pulses is AS — 02 — 01 =0. 



FIG. 5. Transition probability P 2 {t) for the case of two-pulse 
excitation in dependence of time-intervals in units of I/cq. 
The parameters are: A/eo = 0.3, Teo = 10, t/T = 6, lu/eo = 
1. Solid line (a), Aq/eo — 0.19; dashed line (b), Aq/eq — 0.12; 
dotdashed line (c), Aq/eq = 0.25. The phase difference of two 
pulses is AS = 02 — 0i = tt. 


pulses with Gaussian envelopes (HZl) that are shifted on 
the time-interval t between centers of the Gaussians. 

The results for the population of excited state for two 
cases of phase delay between pulses are depicted in the 
Figs. |4] and Figs. [5) As we can see, high-level of the 
population P 2 (t) is realized for two-pulse excitation in de¬ 
pendence from the duration of pulses T and time-interval 
between them t. In Figs. |H if the envelope phase differ¬ 
ence between two pulses is zero, A 6 — 62 — 0i =0, we 
observe that time-evolution of the population strongly 
depends from the amplitude of pulse field and the popu¬ 
lation inversion P 2 {t) — Pi(t) = 2 P 2 {t) — 1 is realized for 
definite time-intervals. If the phase delay is changed to 
TT, i. e. the difference between phases A9 — 62 — 0i = tt, 
the effect of the second pulse reverse the excitation pro¬ 
cess. The typical results are depicted in Figs. [S]for three 
values of the field amplitude. As we see, the curves in¬ 
dicate controllable excitation of qubit state by the first 
pulse and its decay produced by the second pulse. We 
conclude that these results depicted in Figs. [Hand Figs. 


[ 5 ] containing the interference effects between pulses are 
essentially differ from the analogous results for one-pulse 
regime (see, Fig. [T]). 

It is easy to realize from the formula (1551) that depen¬ 
dence on the phase shift in probabilities varies with the 
order of the resonance due to the factor cos(fV(02 — ^i))- 
The order of the resonance reduces the required phase 
difference for observed effects by the order of resonance. 
For example, in the case of second-order resonance the 
effect of probability raising and reducing analogous to the 
result shown on Fig. [5] takes place for A9 = 7r/2, instead 
of TT, as it is for the case of first order resonance. 


E. The case of multi-qubit systems 

Recently, systems composed of several coupled qubits 
(see, e.g., [^) are the subject of much attention. There¬ 
fore, in this subsection we shortly discuss the modifica¬ 
tion of above approach for multi-qubit systems. The con¬ 
crete calculations will be done for two coupled qubits. 
The interaction of two-qubit system with driving held 
have been investigated in theoretical and experimental 
publications (see, e.g., [hOl - f^ l. 

The effective Hamiltonian of multi-qubit system is 
written through the Pauli matrices and tunneling am¬ 
plitudes of the individual qubits 

^ A 

H = (40) 


where 


N 


Ho = - y'5a(t)<5-i“^ + F 


a—1 


2 ^ 

a,b 




(41) 


where time-dependent functions ga(t) are determined by 
Eq.(2) for individual qubits 

ga{t) = ^ea + Aa{t)cOs{uJat +9a) (42) 

and Ja,b are the coupling energies between qubits. 

In this case the system state in the Furry picture 
|4>(t)) = R{t)\^{t)) is calculated through the operator 


R{t) = exp 


i^Aa{t)c 


.(a) _ W ■ 


2 ^ 
a.b 


Ja,b^ 




(43) 


where 

^a(t) = f ga{t')dt', 

JO 


(44) 


in the following form 


|<i)(t)) = Texp 


hmq{t')dt' 


exp(Hm,(t)) (45) 
















through the Hamiltonian of multi-qubit system in Furry 
representation 


where 




( 48 ) 


hmg(t) =(46) 

^ a ' 

and the corresponding Magnus operator Qrngit)- Thus, 
the results of Magnus expantions (12), (13) take also 
place for multi-qubit systems with the Hamiltonian 

^mq (0 ■ 

Below we calculate the Hamiltonian hmq{t) for a two- 
qubit system. In this case 

i? = i?o - i + A2di^)) , (47) 

I 


Thus, the effective Hamiltonian of coupled two qubits in 
the Furry picture is written as 

hmqit) = Ri^), ( 49 ) 


where 


R{t) = exp 


iAi(t)(T^^^ -I- iA2{t)a^^'> — 


(50) 


and J = Ji .2 = J 2 ,i is the strength of the interaction 
between the qubits. 

The transformed Hamiltonian is calculated in the fol¬ 
lowing form 


hmg(^) — Ai 


-fA2 


V 


0 

0 

c2i{Jt+Ki{t)) 

0 


V 

/ 0 
V 0 


^—2i{Jt+K\{t)) 

0 

0 

0 


/ 0 

g2qjt-t-A2(t)) 

0 

0 

0 
0 


g-2i(Jt-|-A2(t)) Q 

0 0 

0 0 

0 g-2i(Jt-A2(t)) 

^^g-2i(Jt-|-AiAi(t)) 

0 
0 


0 

g2j(Ji-Ai(t)) 

0 

0 

0 

0 

g2qjt-A2(t)) 

0 


0 

g2i( Ji—Ai (t)) 

0 


(51) 


and the details of calculations are presented in the Ap¬ 
pendix. This Hamiltonian operates on the four-level sys¬ 
tem of coupled qubits. It should be mentioned that 
this result is obtained in the general form and describes 
two-qubit system interacting with two component driv¬ 
ing field with arbitrary time-dependent amplitudes and 
phases. 


which means independent qubit-pulse interactions for ev¬ 
ery consecutive pulse. The overall effect of pulses in 
this configuration depends on phase changes during the 
pauses. In this case the exponential function reads as 

fc V / 


III. MULTIPLE RESONANT EXCITATIONS 
WITH SEQUENTIAL PULSES 

We focus now on the interaction of the qubit with a 
train of pulses. A more general case is when the external 
field is a sum of pulses with the envelopes Ak (t) and with 
the phases 9k, thus 

1 “ 

9{i) -so + ^Ak{t)cos{ujt-\-9k). (52) 

fc=i 


while the components of the vector p{t) = G{t)/\G{t)\ 
from Eq. (IT^ can be calculated in the following form 


Gxi.t) Ri (-l)^A^cos(A6»fe)^ Jjv dt', 


k^l 

M 


Gyit)^{-l)^Ay2MN9k) [ Jn 

k=i -^0 


f 2Akit') 

V W 


c?t(54) 


G,(t) 0. 


On the whole the probability of qubit state excitation 
reads as 


Below we analyze the effect of a sequence of pulses 
considering the case of not overlapped envelopes Ak{t), 


P2(t) =sin2 (G(t)). 


(55) 
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The obtained results are valid for an arbitrary multi¬ 
resonant interaction, where the order of resonance N is 
indicated by the index of the Bessel function and by 
phase-argument in sine/cosine (see formula (IMll l. These 
results are applicable for an important case naturally re¬ 
alized in practice that is the train consisting of identical 
pulses with equal pulse areas and different phases that 
can be obtained by using various modulators. 


A. Excitation by identical pulses: formation of 
quasienergetic states 

Another interesting configuration that can be realized 
involves the train consisting of identical pulse envelopes 
with the same phases. The trains consisting of big num¬ 
bers of pulses approximately display periodicity in time 
that can be exploited for realization of new regimes of 
qubit dynamics. 

For a sequence of pulses with identical 6k phases G{t) 
can be written as 

^ nt / n A /j.r\ \ 

G(t) = A^ / 

fc=i-^o 


where 


k=l 

This formula is obtained in approximation of not overlap¬ 
ping envelopes, where interactions conditioned by pulses 
are separated by interaction-free time intervals. Thus, 
we get to the equation 

Git) = aJ^Jn dt'. (58) 

To evaluate these quantities further we need to make 
integration of a periodic function Jn ^ ^ ^ with the 

period that equals to time intervals between pulses. For 
train of Gaussian pulses 

OO o 

Ait) = Aq e (59) 

k— — oo 

the period equals to r. 

It is easy to represent the function G(t) as the sum 
of term with linear dependence of time and periodic on 
time term 


Jn ( ) dt' 




M 


Mt) = ^^kit)- 


(57) 


Git) = Aywt + (t)Nit)- (60) 


Here, 4>Nit) is a periodic function with period r defined 
in t G [to Ao + t] as 


Mt) = A 



FIG. 6. Transition probability P 2 (t) for the case of identical 
pulses train in dependence of time-intervals, in units of 1/eo. 
The parameters are: A/eo = 0.3, Teq = 10, uj/eq = 1. Solid 
line (a) describes resonance of the first order, with Aq/eq = 
0.19, r/T = 4; dashed line (b) describes resonance of the 
second order, with Aq/eq = 0.5. 



FIG. 7. Population of excited state P 2 (t) of the qubit under 
sequence of identical pulses showing regular behaviour in the 
case of Ent = tt. Solid line (a) describes resonance of the 
first order, A/eo = 0.3, Teo = 20, r/T = 3 and Aq/eq = 
0.315; dashed line (b) describes resonance of the second order, 
A/eo = 0.4, Teq = 65, r/T = 2.15 and Aq/eq = 0.457. 


and 





(62) 


is the Bessel function averaged on the period. It is easy 
to realize that the last expression can be written through 
the integral of Bessel function for just one of the identical 
pulse envelopes as 



(63) 


The population of the excited states (if the system was 
initially in the lower state, Gi(0) = 1, G2(0) = 0 ) as a 
function of time is then given by 

P 2 it) = sin^ (AyAft -b </>Ar(t)). (64) 

It should be mentioned that the formulae (1591) allows 
to introduce the quasienergetic states and quasienergies 
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for the qubit interacting with time-periodic sequence of 
pulses. Indeed, it is easy to check that the solution of 
Eq. (l20)) for this configuration can be expressed in the 
adiabatic basis 


|±) = ^(|1)±|2)) (65) 

as the quasienergetic states 

|$±) = |±) = |±) ( 66 ) 

with the quasienergies 

E+ = -E- =En = ^J Jn dt. (67) 

The Eq. (I5S1) can be written also as 



Rise arrplitude (in irits cf £b) 


FIG. 8. Quasienergy En depending on the pulse amplitude in 
units of £ 0 . The parameters are: A/eo = 0.3, the pulse width 
Teo = 10, t/T = 4, Ld/eo = 1. Solid line (a) describes the 
first-order resonant configuration; dashed line (b) describes 
the second order resonance. 


where 


|$±)=e^*^-‘C/±(t)|±), 
(7±(t) = 


( 68 ) 


(69) 


is periodic on time in accordance with denotation of 
quasienergetic states 


(7±(t + r) = C/±(t). 


(70) 


It is important also to rewrite quasienergies with the in¬ 
put of qubit energy levels. In this way, taking into ac¬ 
count the relation 




(71) 


we get the states of driven qubit in the following form 



Rise wdth (in nits cf 1/£b) 

FIG. 9. Quasienergy En depending on pulse duration. The 
parameters are: A/eo = 0.4, lj/eo = 1, reo = 40, the ampli¬ 
tude Ao/eo = 0.38. Solid line (a) and dashed line (b) describe 
the cases of first- and second-order resonances, respectively. 


Kit)) = 


p-j J2k ^k{t) cos{bjt+0k) 


| 1 ) 


-h cos{u)t+ 0 k)\^ 2 ) 


(72a) 


Kit)) = e-K+^Ku-{t) 


gieotg-jIDfc Tfe(t)cos(cjt-|-6lfc)|2^^ 


_ g® Yk ^k{t) cos(cjt-l-efc) |2^ 


(72b) 


In these wave functions we added the phase factors in¬ 
volving the sum of qubit energetic levels, that is zero for 
the case of the truncated Hamiltonian ([1]), in which the 
half of the sum of qubit energetic levels has been omitted. 

From the last equations the relation between real en¬ 
ergies and quasienergies can be written as 

Ei^n = ^1 + £')v = ^1 + Ejn, (73a) 

E2 ,n = £2 + Ej^ = £2 — En- (73b) 


Now we analyse time-evolution of the probability of ex¬ 
citation (l64l) . The typical results are depicted in Fig. (|6]) 
for two cases of first-order and second-order resonances. 
As we see the populations for both resonance conhgu- 
rations show complete population inversion for definite 
time-intervals. Such result is cardinally different from 
the case of single-pulse excitation, where the popula¬ 
tion P 2 it) for the second-order resonance is small. It 
can be observed from comparison of Fig. m and Fig. 
©■ The function G{t) = -I- 0Ar(t) is an increasing 

function in time but it grows also periodically due to its 
“linear-l-periodic” structure. Therefore, the dynamics of 
populations Eq. dSl seems to be aperiodic in time as is 
demonstrated in Fig. [SI 

An interesting regime of Rabi oscillations can be 
formed if the following condition for the quasienergy 
takes place 


A'jnt = Ent = 7 rm 


(74) 


Note, that the sum of two quasienergies obeys the rela¬ 
tion E^ -I- E^ = 0, while Ei^n + 7?2,Ar = £i + £ 2 - This 
result is in accordance with the exact result ta king place 
for a two-level atom in a monochromatic field (4l|. 


or in the following form 


A 



Jn 



dt = TTTO, 


(75) 
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where m = 1, 2,.... In this case we observe that G{t+T) = 
A77vt+^Ar(t)+7rm = G{t)+Trm and hence the transition 
probability P 2 {t) = sin^(G(t)) shows periodicity on time- 
intervals. Thus, for the definite values of quasienergy the 
transition probability displays regular, periodic dynam¬ 
ics. The typical results showing clear regular behaviour 
are depicted in Fig. [71 

In the end of this Sec. we analyse the quasienergy in 
dependence from the amplitude of pulses and their du¬ 
ration. Calculations on base of Eq. (1571) lead to large 
values of the quasienergy in comparison with the anal¬ 
ogous results for two-level atomic system. As we see in 
Fig. [8] and Fig. [9] the quasienergy increases almost lin¬ 
early with increasing the amplitude of pulses as well as 
with the duration of pulses for both resonant configura¬ 
tions. 


This approach has been formulated in general case for 
analysing the Rabi model with time-dependent interac¬ 
tion Hamiltonians along the z axis for multi-resonant in¬ 
teractions. Nevertheless, we demonstrate that the ap¬ 
proach proposed can be used for pulse-driven multi-qubit 
systems as well as for interaction of qubits beyond RWA. 
Thus, an interesting question to address would be appli¬ 
cation of the general expressions obtained in this paper 
for qubit operations beyond RWA and for elaboration of 
quantum gates with several qubits. 


ACKNOWLEDGMENTS 

We acknowledge support from the Armenian State 
Committee of Science, the Project No.I3-IC03I. 


IV. CONCLUSION 


Appendix 


In conclusion, we have investigated the dynamics and 
manipulations of qubit by its multiple resonant interac¬ 
tions with external pulses. Single qubit quantum op¬ 
erations have been proposed and analyzed for various 
regimes of state excitation including: single-pulse excita¬ 
tion, two-pulse excitation with phase shift between pulse 
envelopes being controlling parameter and for excitation 
with identical sequential pulses. This investigation have 
been performed in the approach, based on Bloch vec¬ 
tor’s treatment and the Magnus quantum electrodynam¬ 
ical formalism in Furry picture that can be applied for 
various types of interactions. 

This work has showed the possibility to produce a large 
variety of qubit operations for time-domain in depen¬ 
dence of the amplitude of driven field, the parameters of 
pulses, the tunneling rate and the controlling phase fac¬ 
tors. Considering two resonant configurations, i.e. first- 
and second-order resonant configurations, we have shown 
that the complete population inversion can be achieved 
through tailoring the pulse train for both configurations 
at definite time intervals. 

Particular emphasis has been laid on interaction of 
qubit with train of identical pulses. We have shown that 
quasienergies and quasienergetic states of combined sys¬ 
tem ” qubit -|- pulse train” are forming for all multiple res¬ 
onant regimes. The quasienergy increases almost linearly 
with increasing the amplitude of pulses as well as with 
the duration of pulses. Moreover, we have demonstrated 
that for all multiple resonant regimes the state popula¬ 
tion exhibits time-dependent aperiodic oscillations, but 
also displays periodically regular oscillations for the def¬ 
inite values of the quasienergy, E]\[T = nm. 


To calculate the Hamiltonian (I49p we use the following 
well-known formulas 


e*“'^"'(Ta,e cos(2Q!)(Ta, — sin(2Q!)(Ty, (A.la) 

= sin(2Q!)(Tx-b cos(2a)(Ty, (A.lb) 

that allows us to performe unitary Furry transfomations. 
In this way, for the Pauli matrice one of the qubit we 
obtain the following operator equation 


= cos(2Ai(f)) cos{2Jta^p) — sin(2Jti7^^^)^ 

— sin(2Ai(t)) sin(2JttT^^^) -b cos(2 . 

(A.2) 


Then, this equation is transfomed in the matrice form as 
the following 


R-\t)dlR{t) 


0 

g2i g2iAi (t) 


g-2jJta(’')g-2jAi(t) 


/ 0 0 

0 0 

g2i(Ji-|-Ai(i)) Q 

V 0 g-2i(Jt-Ai(t)) 


g-2i(Jt+Ai(t)) 

0 

0 

0 


e 


0 \ 

2j(Jt-Ai(i)) 


0 

0 


(A.3) 

) 


The expression R R(t) is evaluated in a similar 

way. On the whole, these results lead to the formula m- 
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